Quasi-one-dimensional disordered systems: fluctuations, transport and interplay 
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In a one dimensional lattice thermal fluctuations destroy the long-range order making particles 
of the lattice move on a scale much larger than the lattice spacing. We discuss the assumption that 
this motion may be responsible for the transport of localized electrons in a system of weakly coupled 
chains. The model with diffusing localization sites gives a temperature-independent mobility with 
a crossover to an activated dependence at high temperature. This prediction is consistent with and 
might account for experimental results on discotic liquid crystals and certain biopolymers. 
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I. INTRODUCTION 

In one and two dimensional (D) systems positional cor- 
relations diverge with the system size due to thermal 
fluctuations. Such loss of long-range order, known as 
Landau-Peierls instability, implies the divergence of the 
mean square displacement (MSD) of system's structural 
units. For harmonic lattices this result may be illustrated 
by direct calculations. Let qi be a displacement of the z-th 
atoms in a harmonic chain from its equilibrium position 
Xi ~ ia where a is the lattice spacing. Assuming that 
the atoms are distributed as in thermal equilibrium, one 
can show that the atomic MSD is proportional to the 
atom's distance from the chain's end, 



k 2 > = 



(1) 



Here u>q = yjk/m, k is the harmonic force constant, 
and m is the mass of an atom. In a sufficiently long 
chain the atomic displacement from equilibrium position 
(g?) 1 / 2 may be significantly longer than the equilibrium 
lattice spacing a. For instance, for the force constant 
k ~ 1 N/m, temperature T ~ 10 2 K, and i = 10 4 the 
equation (1) predict the displacement of order 100 A. 
In the limit of the infinite chain the MSD diverges and 
atomic motion is unbounded. The loss of long-range or- 
der due to Landau-Peierls instability has been observed 
in X-ray and neutron scattering experiments in many 
quasi- ID and quasi-2D systems such as liquid crystals 
and membranes. 

The divergence of the MSD does not mean, of course, 
that the chain is unstable. While the long range order is 
lost, the short order is preserved. Indeed, the result for 
the relative displacement of two atoms reads 



((qi - q i+3 ) 2 ) 
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(2) 



According to this equation the distance between two 
adjacent atoms (j — 1) does not depend on the 
chain's length and deviates from a by a value of order 
(fcflT '/mwj) 1 ' 2 , which is normally small. 



In recent paper [lj it was suggested that since individ- 
ual atoms in a ID lattice are to some extent delocalized, 
they may serve as temporary vehicles for localized elec- 
trons. Consider a system of parallel chains separated by 
the distance b which is larger than the lattice spacing a. 
As will be discussed in the following sections, on a long 
time scale (u>ot ^S> 1) a tagged atom in a long isolated 
harmonic chain behaves as a Brownian particle with the 
diffusion constant 



D„ 



kT 
2mu>o ' 



(3) 



which increases linearly with temperature. On the other 
hand, the diffusion constant Dh for electronic hopping in 
a static disordered chain depends on temperature expo- 
nentially and may be smaller than D a for sufficiently low 
T. For instance, for the nearest-neighbor hopping over 
uncorrelated sites with the Gaussian energy distribution 

ff (e) = C exphe 2 /2e 2 ], (4) 

the hopping diffusion constant has the form 



D h = — v exp 



-a/L- (e /k B TY 



(5) 



where L is the localization length of the carrier wave 
function, and v is the "attempt frequency" . For the typ- 
ical parameter set k ~ 10 N/m, eo ~ 0.1 eV, a ~ 10 A, 
v/ojQ ~ 1, and exp(— a/L) ~ 10~ 3 , one finds that 
Dh <C D a when eo/ksT > 1. This estimation suggests 
what we call the hitchhiking mechanism of electronic 
transport: for sufficiently low temperature the hopping 
mechanism is responsible only for the transport perpen- 
dicular to the chains, while along the chains electrons are 
transported predominantly by mobile localization sites. 

To estimate the diffusion constant Dh for the lateral 
transport along the chains due to the hitchhiking mecha- 
nism let us assume that the localization is strong and each 
atom is associated with a mobile localization site whose 
MSD as a function of time (Aq 2 (t)) = ([q(t) - q(0)} 2 ) 
is given. Suppose also that the transition rates for the 
inter-chain hopping do not strongly fluctuate around a 
typical value W± . Then the lateral motion of an electron 
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can be considered as a ID random walk with the the time 
step r ~ l/W± and the length step I ~ (Ag 2 (r)) 1 / 2 . The 
corresponding diffusion constant is 1 2 /2t, which gives 



(6) 



If the atomic MSD is diffusive (q 2 (t)) ~ 2D a t, as in an 
isolated chain, then £)ii does not depend on W± and co- 
incides with the atomic diffusion coefficient 



D\\ = D a 



k B T 



(7) 



2mwo 

The corresponding mobility is temperature independent 



eD 



fi 



2k ' 



(8) 



For uiq — 10 13 s _1 and k — 10 N/m this equation gives 
p ~ cm 2 (V / s)" 1 , which is consistent with experi- 
mental values for columnar liquid crystals and certain 
conjugated polymers. 

Temperature independent mobility was observed in 
many low-dimensional soft matter systems, in particular 
in discotic liquid crystals [3| and DNA J4j . It is usually 
explained in terms of the polaron model [3J] or the dynam- 
ical disorder models [5[ . The model of hitchhiking trans- 
port gives an alternative and very simple explanation, 
which does not involve any adjustable parameters. The 
model also predicts a crossover to an activated temper- 
ature dependence for sufficiently high temperature when 
the inequality Dh <C D a is no longer valid, and the con- 
ventional hopping mechanism begins to dominate. Such 
crossover was reported for charge transport in DNA Q ■ 

The underlying issue of the model is the one of delocal- 
ized single-particle dynamics in low-dimensional lattices. 
The purpose of this paper is to illustrate this issue for two 
simplest approximations, one of non-interacting chains, 
and the other of dissipative chains, subjected to the noise 
and friction forces. 



II. LACK OF LONG-RANGE ORDER 

Although derealization of an atom in low-dimensional 
lattices is a well-known result, it might appear counter- 
intuitive and is often misinterpreted. For instance, in Q 
the the divergence of the atomic MSD in a ID chain 
is considered as an artifact which arises from the zero- 
frequency mode in the chain with periodic boundary con- 
ditions. Actually, derealization is the general property 
of low-dimensional lattices and does not depend on the 
type of boundary conditions. The question was compre- 
hensively studied by Montroll [7J] , but the generality of his 
approach makes it rather complicated. In this section we 
give a simple derivation of Eq. (1) for a harmonic chain 
and generalize it for the quantum case. 

Consider a harmonic chain of N + 2 atoms with fixed 
ends. Labeling atoms by index i = 0,1, ...N + 1 one can 



write the Hamiltonian in the form 



H 



N 9 N+l 



(9) 



assuming that displacements for the terminal atoms are 
zero, qo = qN+i = 0. The Hamiltonian can be diagonal- 
ized by means of a normal mode transformation, 



1 N 



Pi 



3=1 



N 

3=1 



with normalized eigenvectors 
Al] = ( ATI 



sin 



N + l 



(10) 



(11) 



which satisfy the orthogonality condition J3j=i ^-ij^ik = 
5jk- In terms of normal coordinates the Hamiltonian as- 
sumes the form, 



3 = 1 

where the normal mode frequencies are 

U) a = 2LUn sin — 

3 |_2(A + 1)_ 

Then the atomic MSD can be written as 

33' 



(12) 



(13) 



(14) 



where the average is taken with respect to the canonical 
distribution p — Z~ 1 e~ H / kBT . Since 

(QiQi>) = %' (PjPi>) = S jS >k B T, (15) 

U 3 

the equation (TT4"]) takes the form 

N A 2 



knT v— -\ A 



E 

3 = 1 



(16) 



Inserting the expressions for Aij and u>j , and introducing 
the new variable 



2(N+1) 



one obtains 



k B T 



N 



ft = 



4mw 2 N - 



sin 2 (2i6j) 



3 = 1 



Sill 



(17) 



(18) 



Since A0j = 9 3+1 - 6j = n/2(N + 1), the sum in the 
above expression can be converted in the limit N — > oo 
into the integral as follows 



muj 2 ir J sin 2 ( 



(19) 
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which eventually gives 



k B T . k B T . 



(20) 



The equation (2) can be derived in a similar way. 

In fact, the result (|2T)|) can be obtain without the nor- 
mal mode transformation, but using instead new coor- 
dinates Si = qt — q-i-i and noticing that qi — Y^k=i&k- 
This formula is the manifestation of cumulative nature 
of deviation from equilibrium in ID systems. Another 
comment concerns the equations (|16[) . the summation of 
which gives 



£k 2 > =— £< 

* — * in — ^ 



(21) 



This relation does not involve Aij and is believed to be 
a general result [8J]. 

The shortest way to generalized the result (|20|) for the 
quantum case is to use the quantum mechanical formula 
for the average energy of the oscillator, corresponding to 
a j-th mode 



/ fadj \ 



^•> = ¥ coth U^> (22) 

Since the average potential and kinetic energies are equal, 

f($) = 2< p /> = 2 {Ej) = Y coth [r 



2k B T 



one obtains 



(QjQf) = 5ji> - — coth ' 



2u!i 



\2k B T 



Substitution of this result into (Q3J) gives 

Hluq sin 8 



, 9l h 1 sin 2 (2i0) 

(» > = / r } ~ S - L 

muio ir J sm a 



coth 



k R T 



(23) 



d0. (24) 



In the high temperature limit, hujQ/k B T -C 1, one can 
use the approximation coth(x) ~ 1/a;, which leads to the 
classical result (|2U|) . 

For the ultimate quantum case T = 0, coth(x) goes to 
one, and the MSD takes the form 



mujQ 7r J sm v 
where the integral increases with i logarithmically 



5 sin 2 (2£0) 



1 



dO ~ -InVi 
sin 2 



(25) 



(26) 



Thus the MSD due to quantum zero-point fluctuations 
reads as follows, 



2irmuJo 



In z. 



(27) 



For luq ~ 10 12 s 1 , to ~ 10 27 /eg (proton), and lni ~ 1, 
the above equation gives (q 2 ) 1 ^ 2 ~ lA. 



III. DYNAMICS 

With the delocalized character of atomic motion in a 
long chain established, let us consider the question about 
the dynamics of this motion. It can be conveniently de- 
scribed in terms of a velocity correlation function 



(28) 



For instance, integrating C(t±,t2) one obtains the MSD 
(Aq 2 {t)) = ([q(t) - q{0)} 2 ) of an atom: 



(Aq 2 (t)) = [ dt 2 [ dhCituh). 
Jo Jo 



(29) 



Using stationarity of the process v(t), C(ti,t 2 ) = C{t% 
t\) and integrating by parts, one gets 



(Aq 2 (t)) =2 / dr(t-r)C(r). 



(30) 



Next, one can show [l[ that the Laplace- Fourier trans- 
form of the velocity correlation function C(w) = 
/ °° dt e~ luJt C(t) determines the dynamical mobility /Lt(u>) 
of a charged atom, 



k R T 



CV). 



For an isolated chain, C(t) is a Bessel function 



C(t) 



knT 



J {2uj t). 



(31) 



(32) 



This result can be obtained using the normal mode trans- 
formation of the previous section. Indeed, since 



1 N 



(33) 



and 



Pj(t) = Pj{0) cosuijt - ujjQjiO) siaujt, (34) 
one obtains for the correlation Ci(t) = (vi(0)vi(t)) 

1 - 

City) = ~J2 ^4^(0)i\(0)) cosw,-*. (35) 

Assuming that initial distribution of coordinates and mo- 
menta is canonical, one gets (Pj(O)Pfc(O)) = Sjkk B T. 
Then 



Ci(t) 



kuT 



N 



£ A 2 j cos Ujt. 



(36) 



Substituting the explicit expressions for and u>j and 
converting the sum into an integral one obtains 

k T 4 f 7 ^/ 2 

CAt) = — / d9 sin(2i0) cos(2w sin0), (37) 

TO 7T J 
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Ci{t) = — { Jo(2uj t) - J 4l (2w i)}. (38) 

For large i the term with Ja may be neglected , and one 
recovers the result (j3"2")l . 

Using ([30]) and ((32j) one obtains for the MSD 

(Aq 2 (t)) = 2D a t - 2D a tJ 1 {2u Q t), (39) 

were the diffusion coefficient is 

k B T 



D„ 



dtC{t) 



(40) 



For long time, u>at ^> 1, the first term in (|39|) dominates, 
so that the atomic motion is diffusive, (Aq 2 (t)) ~ 2D a t. 
This result was first discussed by Rubin [9j. fn the same 
limit the mobility is purely real and does not depend on 
frequency, fi(u>) w e/2u)am. 



IV. DISSIPATIVE CHAINS 

In previous sections the interaction of chains has been 
neglected which is very likely to be an over-simplification. 
The problem of dynamics of interacting chains may be 
considered as a generalization of the Frenkel-Kontorova 
model about a chain in an external spatially periodic po- 
tential field [ioj . In quasi-lD systems this periodic poten- 
tial is created by adjacent chains and is not static, which 
makes the problem very difficult [ll[. It was suggested 
in [l[ that some insight can be achieved by modelling the 
chains interaction using the Langevin approach. Namely, 
one may assume that the force exerted on an atom by 
adjacent chains can be written as the sum of a regular 
dissipative force linear in the atom's velocity, —7%, and 
a fluctuating term £j (t) . With this assumption the equa- 
tion of motion of the atoms takes the form 

mq\(t) = + q i+1 - 2qi) - jq^t) + &(t). (41) 

Let us treat the fluctuating term as a zero centered 
white noise which is not correlated for different atoms 
and related to the friction constant 7 through the con- 
ventional fluctuation-dissipation relation, 

(£(*)> =0, {^mk(t))=2k B TjS ik S(t). (42) 

Such approach is common in polymer physics to describe 
the polymer-solvent interaction. In that case the fric- 
tion constant 7 can be expressed in terms of the solvent 
viscosity. In our model, the dissipation and fluctuating 
terms describes interaction with other chains, and 7 is 
an adjustable parameter. 

Let us find the velocity correlation function of an atom 
for the model of dissipative chains. Referring again to 
the normal mode transformation (jfOp . the inverse has 
the form 



N N 



(43) 



Then from the above equation of motion one finds the 
equation for Qj(t) 



Qj(t) - -utQjit) - XQj(t) + vM (44) 



Here 



(45) 



has a meaning of the random force for a normal mode 
Qj, and A = 7/771 is the inverse velocity relaxation time 
for an atom. Using the method of Laplace transform the 
solution of (|4"4"|) can be written as 

Qj(t) = aj(t)Qj(Q) + &i(t)Qj(0) + / drb^t - t)^(t). 



(46) 

The explicit form of the functions cij (t) and bj (t) depends 
on the sign of the difference u)j — A/2. We consider here 
only the case of an over-damped chain when A > 2u>o 
(and therefore A > 2u>j for any normal mode j). In this 
case 

a,j(t) = e^T* smhfij/; + cosh^ij , (47) 



bj(t) = —e 2 f sinhOjf, 



(48) 



where flj = yA 2 /4 — uo 2 . Then for the velocity correla- 
tion function C(t) = (qi(t)qi(0)) one gets 

k B T 



C(t) 



(49) 



For the very strong dumping A 3> loq, bj(t) can be ap- 
proximated as 

bj(t) = exp(-Xt) - (ujj/X) 2 exp (-Wjt/\) . (50) 

Transforming the sum (|49p into an integral one obtains 
the velocity correlation function in the form 



C(t) 



k B T _ xt k B T d 

e -I — 

m mX dt 



{e~ at I (at)} (51) 



where Io{x) is the modified Bessel function, and a = 
2u 2 /X. Note that a < A. 

Substituting the above expression for C(t) into (|50"|) 
one obtains 

(Aq 2 (t)) = ^L t e at {l (at)+I 1 (at)}. (52) 

For long time, t 3> a -1 , this expression gives sub- 
diffusive behavior 



(Aq 2 (t)) ~ 2FVt 



(53) 
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with the mobility factor 



Sub-diffusive motion with the MSD growing as yr is 
characteristic feature for any form of the single file dif- 
fusion when overdamped Brownian particles constrained 
to move in one dimension and are not allowed to pass 
each other [T^]. Its consequence for the hitchhiking 
model is that the diffusion constant depends on the inter- 
chain transition rate W±, as follows from equation (JQ, 
D\\ — Fy/W±. Note however, that this result holds only 
for the transport on the time scale much longer than 
1/Wj_; it does not apply for short chains when inter-chain 
transitions are negligible. 

Another consequence of sub-diffusive transport is that, 
in contrast to the approximation of non-interacting 
chains, the dynamical mobility essentially depends on fre- 
quency. Using (f3"Tj) one can find that for oj <C a both real 
and imaginary parts increases with u as oj 1 / 2 , namely 

fih,,) = JL l-^-n-i). (55) 

Power frequency dependence fJ-(u>) ~ oj s with < s < 1 
is typical for many disordered systems. Note however, 
that for frequency lower than the inter-chain transition 
rate, one has to take into account inter-chain hopping of 
carriers. On this time scale the carrier diffusion is normal 
(Aq 2 (t)) ~ t, and is almost frequently independent. 



V. CONCLUDING REMARKS 

This paper promotes the idea that in quasi-lD sys- 
tems long-range structural fluctuations may carry local- 
ized electrons over a considerable distance and that for 
sufficiently low temperature this may be the dominat- 
ing mechanism of charge transport. Most of the results 
are obtained under very idealized assumptions and can 
hardly be regarded as anything but toy-model calcula- 
tions. On the other hand, qualitative predictions of the 
model seem quite general. For instance, while the charac- 
ter of time dependence for the atomic MSD may depend 
on many factors, the linear dependence of the MSD on 
temperature is a general property, which is responsible 
for a temperature independent hitchhiking mobility. 

The assumption of one-dimensional dynamics of the 
chains is also not essential. In this paper we assumed that 
the chains are parallel and form strongly anisotropic 2D 
or 3D crystal. One may argue that the model of hitch- 
hiking electronic transport may be relevant to polymer 
systems, where monomers move sub-diffusively in three 
dimensions 13]. Another generalization is the case when 
the electronic localization length is much larger than the 
lattice spacing, L 3> a. This problem is related to the 
dynamics of a large cluster [l4| rather than of a single 
particle in a chain. 
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